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Abstract 

We propose the multi-state complete-active-space second-order perturbation theory spin-orbit 
method (MS-CASPT2-S0) for electronic structure calculations. It is a two-step spin-orbit coupling 
method that does not make use of energy shifts and that intrinsically guarantees the correct 
characters of the small space wave functions that are used to calculate the spin-orbit couplings, in 
contrast with previous two-step methods. 
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I. INTRODUCTION 



In electronic structure two-step spin-orbit coupling methods, dynamic correlation is han- 
dled in the first step, using the spin-free part of the Hamiltonian and a large configurational 
space in variational or perturbational schemes. Then, spin-orbit coupling is handled in the 
second step, using an effective Hamiltonian and a small configurational space in spin-orbit 
configuration interaction (CI) calculations. In these methods, the effective Hamiltonian con- 
tains explicit energy shifts, which are a mean to transfer dynamic correlation effects from 
the first step to the second step in a simple and effective manner. 

It has been found that the energy shifts of the spin-orbit free levels, which are driven 
by their energy order within each irreducible representation, can lead to anomalous results 
when avoided crossings exist with significant change of character of the wave functions at 
each side, which take place at different nuclear positions in the large and in the small 
electronic configurational spaces. In these cases, the shifts must be assigned according to 
the characters of the wave functions.^ This usually implies analyses of wave functions in 
both configurational spaces. 

The ultimate reason behind these problems, which are present in the available two-step 
methods, is the different nature of the wave functions of the spin-free states in the large 
and in the small configurational spaces, so that, even when the avoided crossings do not 
exist or when they take place at the same nuclear positions in the large and in the small 
spaces, such a different nature makes the spin-orbit couplings calculated in the small space 
not as accurate (meaning as close to the spin-orbit couplings calculated in the large space) 
as desired. 

Here, we propose the multi-state complete-active-space second-order perturbation theory 
spin-orbit method (MS-CASPT2-S0). It is a two-step method that does not make use of 
energy shifts and that guarantees by construction the correct characters of the small space 
wave functions that are used to calculate the spin-orbit couplings. 
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II. THE MS-CASPT2-SO METHOD 



Let us assume we have a many electron system with a Hamiltonian H which is made of 
the addition of a spin- free contribution, H^^ , and a spin-orbit couphng contribution, H^^: 

H = H'^ + H^"" . (1) 

In the spin-orbit free MS-CASPT2-S0 method, the procedure is initially the same as 
the MS-CASPT2 procedure:^'® Several state-average complete-active-space self-consistent- 
field SA-CASSCF (or CASCI) states are calculated, which define a reference con- 
figurational space called the V space. Let us collect them in the row vector 
^CAs ^ Q ^CAS^^ I ^CAS^^ J where p is the total number of SA-CASSCF states. 

These wave functions can be classified according to their values of spin quantum numbers 
and symmetry group irreducible representations and subspecies, SMsT^, but we will omit 
these labels here for simplicity. 

In spin-orbit free MS-CASPT2 calculations, the SA-CASSCF wave functions are used as 
a basis to calculate the matrix of a spin-free second order effective Hamiltonian, 
which is defined in Eq. 30 of Ref. 5 and depends only on the spin-free part of the Hamiltonian, 
H^^. This matrix, which is H^nd^^'^^^ = ^^"^^"^ H2nd^^^ '^^^ is diagonalized in order to 
compute the MS-CASPT2 energies, as the eigenvalues, and the modified SA-CASSCF (or 
CASCI) states, as the eigenfunctions: 

Hl^ff^CAS^^^^MS2^ (2) 

where E^^^ is a diagonal matrix with the MS-CASPT2 energies Ef^^, . . . , E^^^^, as 

the diagonal elements and L[ is a unitary transformation of the original SA-CASSCF wave 
functions that preserves the SMsT'-f values, 

^CAS' ^-^CAS^^ 

Obviously, the modified SA-CASSCF wave functions ^'-^^^ also span the V space. What is 
important is that they are the appropriate zeroth-order basis for a second-order perturbation 
theory treatment of the dynamic correlation that leads to the MS-CASPT2 energies^ and 
they have the appropriate characters in correspondence with these energies. 

In the MS-CASPT2 spin-orbit calculations proposed here, we can follow two alternative 
procedures that lead to the same result. Both of them are based on the use of the spin- 



dependent effective Hamiltonian that results from the addition of the spin-orbit couphng 
operator to the spin-free effective Hamiltonian of the MS-CASPT2 method, 

H^^^ = Hlf'' + ■ (4) 

In the first procedure, which is a formal two-step procedure, the regular spin-orbit free 
MS-CASPT2 calculation is completed and the modified SA-CASSCF wave functions '^^^ 
are used as a basis for the matrix representation of H^^^. The resulting matrix, 

^eff,CAS' ^ ^CAS'^ jjeff ^CAS' ^ ^MS2 _^ ^SO,CAS' ^ 

with H^o,CAS' ^ ^CAS't fyso^cAS'^ couples the modified SA-CASSCF states via spin-orbit 
coupling. (It couples different SMsTj blocks and it can be factorized according to double 
group irreducible representations.) Its diagonalization gives the final energies and wave 
functions: 

^eff,CAS' jjSO' ^ jjSO' ^MS2-30 ^ ^g-j 

where E_ ~ is a diagonal matrix with the MS-CASPT2-S0 target energies 
^MS2-so ^^MS2-so ^ _ _ _ ^^Msi-so ^^le diagonal elements and Uf°' is a unitary trans- 
formation of the modified SA-CASSCF wave functions that couples the SMsT'^ values and 
gives the target spin-orbit wave functions, 

^M52-SO ^ -qiCAS' jjSO' 

Alternatively, in the second procedure, which is a formal one-step procedure, the original 
SA-CASSCF wave functions '^'^^ are used as the basis for the matrix representation of 
H'^^f . In order to do this, the regular spin-orbit free MS-CASPT2 calculation does not need 
to be completed, but only the computation of the H.2nd^^''"^^ matrix used in Eq. 2, plus 
the addition of the matrix of H^'^ in this basis (H^o,cas ^ ^CASf ^so ^cASy 

^eff,CAS ^ ^CASj jjeff ^CAS ^ jjSF^eff,CAS _^ jjSO,CAS ^g-j 

Its diagonalization gives the same target energies and wave functions as the first procedure, 

^eff,CAS jjSO ^ ^SO ^MS2-SO ^ j^g-j 

with 

^MS2~SO ^ \^CAS jjSO ^20) 



The present method is closely related with the restricted active space state interaction 
approach with spin-orbit coupling of Ref. 3, SO-RASSI. The latter, when it is used in a 
SA-CASSCF/MS-CASPT2 context, corresponds to diagonalizing the mixed effective Hamil- 
tonian matrix E_^^^'^ + }jf'0,CAS ^ with the spin-orbit free part of Eq. 5 and the spin-orbit 
coupling part of Eq. 8. The results of both approaches are expected to be similar when the 
CAS and the CAS' wave functions (Eq. 3) are also similar, which is a common case. Differ- 
ences should show up when the two sets of wave functions are not so similar in a one-to-one 
basis, for instance when the dynamic correlation switches the order of states. This is the 
basic advantage of the present method. However, we must note that the SO-RASSl method 
can be used together with general single-state and state-average RASSCF and CASSCF plus 
RASPT2 and single-state and multi-state CASPT2 spin-orbit free frameworks, whereas the 
present one can only be used in the spin-orbit free framework of SA-CASSCF plus MS- 
CASPT2 calculations. 

Let us now justify why H'^-^^ (Eq. 4) is the effective Hamiltonian of choice in the MS- 
CASPT2-S0 method. For this purpose, we recall that the basic idea of two-step methods is 
to use a spin-orbit effective Hamiltonian made of a spin-free effective Hamiltonian (usually 
ff'^^ + i^'^^*-^*) plus the spin-orbit coupling operator H^'^, and to choose the spin- free effective 
Hamiltonian by imposing the requirement that, when used in the small space V of the second 
step, it has the same eigenvalues that H^^ has in the large space Q of the first step.^ In 
the particular case in which the first step is a MS-CASPT2 calculation and the small space 
of the second step is defined by the SA-CASSCF wave functions, H2nit^^ is a spin-free 
effective Hamiltonian that fulfills such a condition. In consequence, W^^^ (Eq. 4) is the 
proper spin-orbit effective Hamiltonian. 

III. CONCLUSION 

A two-step spin-orbit coupling method for multi-state complete-active-space second-order 
perturbation theory calculations MS-CASPT2 is proposed which does not make use of energy 
shifts. It intrinsically guarantees the correct characters of the small space wave functions 
used to calculate the spin-orbit couplings, in contrast with previous two-step spin-orbit 
coupling methods, where it has to be checked externally. 



5 



Acknowledgments 



This work was partly supported by a grant from Ministerio de Ciencia e Innovacion, Spain 
(Direccion General de Programas y Transferencia de Conocimiento MAT2008-05379/MAT). 



^ R. Llusar, M. Casarrubios, Z. Barandiaran, and L. Seijo, J. Chem. Phys. 105, 5321 (1996). 

2 V. Vallet, L. Maron, C. Teichteil, and J.-P. Flament , J. Chem. Phys. 113, 1391 (2000). 

3 P. A. Malmqvist, B. O. Roos, and B. Schimmelpfennig, Chem. Phys. Lett. 357, 230 (2002). 
^ G. Sanchez-Sanz, Z. Barandiaran, and L. Seijo, Chem. Phys. Lett. 498, 226 (2010). 

^ J. Finley, P.-A. Malmqvist, B. O. Roos and L. Serrano- Andres, Chem. Phys. Lett. 288, 299 
(1998). 

6 A. Zaitsevsisii and J. P. Mahieu, Chem. Phys. Lett. 233, 597 (1995). 



6 



